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Abstract 

Four-dimensional heavy-fermion QED is studied in light-cone coordinates with 
(anti-)periodic field boundary conditions. We carry out a consistent light-cone canon- 
ical quantization of this model using the Dirac algorithm for a system with first- and 
second-class constraints. To examine the role of the zero modes, we consider the quan- 
tization procedure in the zero-mode and the non-zero-mode sectors separately. In both 
sectors we obtain the physical variables and their canonical commutation relations. 
The physical Hamiltonian is constructed via a step-by-step exclusion of the unphysical 
degrees of freedom. An example using this Hamiltonian in which the zero modes play 
a role is the verification of the correct Coulomb potential between two heavy fermions. 
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1 Introduction 

There has been recent progress [1-4] towards a framework for describing decay pro- 
cesses involving heavy quarks. Its basis has been the study of QCD in which some of the 
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quark masses are taken to infinity. In this limit, the quark spin degrees of freedom decou- 
ple from each other and the couplings of the heavy quarks to the gluon degrees of freedom 
become mass independent and are described by a Wilson line. 

The remarkable consequence of this independence is that all matrix elements of vector and 
axial vector currents bilinear in either or both of a pair of heavy quarks taken between either 
vector or pseudoscalar mesons in initial and final states at arbitrary momentum transfer 
are determined in terms of a single normalized function. The question arises as to how one 
calculates this function. 

One of the more promising approaches to problems in QCD is light-cone quantization 
1^, ^. Light-cone quantization has turned out to be a useful tool for the perturbative 
treatment of field theories 0, In its extension [9-13] to the nonperturbative domain, 
we have come to realize that careful attention must be paid to the non-trivial vacuum 
structure of light-cone quantum field theory. For instance the light-cone vacuum in the 
massless Schwinger model can be understood only by a careful study of the zero-modes of 
the constraints imposed by the light-cone frame [14-16]. Indeed it has been conjectured that 
the dynamics of zero modes in QCD in light-cone quantization provide the mechanism for 
confinement ||, |^ . 

In the present paper we consider a four- dimensional heavy-fermion QED, as an initial step 
towards the study of heavy-quark QCD in light-cone quantization. We carry out a consistent 
canonical quantization of this gauge invariant model in a light-cone domain restricted in its 
"spatial" directions. It is well known that in such a restricted region one has problems 
with zero modes [17,15,16]. The canonical quantization of the massless Schwinger model 
on a light-like hyperplane, taking into account zero-mode contributions, was considered in 
15, 1^. Neglecting such contributions. Tang et al. |18| carried out the discretized light-cone 



quantization of four-dimensional QED and Mustaki [|T9| developed the canonical quantization 
of two-dimensional QED on the null plane. We note that these questions have also been 
investigated by first quantizing and then taking the heavy-quark limit, in two-dimensional 

QCD g. 

In studying the quantization of gauge field theories one is confronted by first-class con- 
straints and their corresponding gauge conditions. A consistent canonical quantization for- 
malism for such problems was proposed long ago by Dirac [^] and Bergmann and its 
generalization to fermionic (Grassmann-odd) constraints by Casalbuoni (There are in 
fact books on this subject ||2^, We will not discuss here other approaches involving path 
integration, BRST or BFV methods.) In the Dirac approach some problems arise. They 
involve the determination of the dynamical (physical) variables and the construction of the 
physical Hamiltonian. One needs to prove, also, that in the physical sector the S'-matrix 
does not depend on the form of the gauge conditions. 

The specific gauge theory we address in this paper is in terms of light-cone variables 
where, as we shall see, the quantization comes with some important constraints involv- 
ing zero-frequency variables which require special attention. To examine the role of the 
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zero modes explicitly we consider the quantization procedure in the separated zeroless-mode 
and zero-mode sectors. In such sectors we choose special gauge conditions and obtain the 
physical variables and their canonical commutation relations. The physical Hamiltonian is 
constructed by systematic elimination of the nondynamical variables. 

In Sections 2-5 the canonical quantization of heavy fermion QED is addressed by way of 
the Dirac-Bergmann algorithm. In Section 6 we calculate the interaction potential between a 
heavy fermion and a heavy antifermion using old-fashioned perturbation theory. Conclusions 
and a brief discussion follow in Section 7. 



2 Constrained Dynamics of Heavy Fermion QED 

We begin with the Lagrangian of heavy-fermion QED obtained by a generalized 
Foldy-Wouthuysen transformation that removes the terms mixing fermion and anti-fermion 
fields in the action. We have 

£ = i^I/IU^'D^'^ - il/M^ - ^F^^F^"", (2.1) 

where the Minkowski metric is diag^f^j^ = (1,-1,— 1,-1), 

D^ = df, + leA^ (2.2) 
is the covariant derivative, 

JjKJ^ = UU^, (2.3) 

and is the given 4- velocity of the heavy fermion satisfying the condition [/^ = 1. The 
field strength tensor is 

F^u = d^A, - d,A^, (2.4) 

and we use the system of units where h = c = 1. The heavy fermion limit means that M[/^ 
is greater than any other momentum in the problem under consideration. 

The light-cone coordinates in four- dimension are = {x~^,x~ ,x^), j = 1,2, where 

x^ = ^{x^±x^) (2.5) 

The variable x^ plays the role of the "time" variable. In terms of such coordinates the 
Lagrangian C becomes 

C = i^JJJ{U+d^ + U^d+ - U.jd^)^> - M'^m + ]^Fl_ + - ^FjkFjk 

-e^l/I{U+A_ + f/_A+ - UjAj)-^, (2.6) 
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where 



Fjk — djAk — dkAj, 
F+_ = d+A^ - 
F±j = d^Aj-djA±, 

d± = d/dx^ (2.7) 

The fermion part of the Lagrangian ( |2.1| ) differs from ordinary QED by the change 
7^ — > JJJU'^. Let us consider the equations of motion for the fermion field They are 

iJjJU^'D^m = (2.8) 

As in ordinary QED, let us introduce, in light-cone coordinates, two "chiral-like" fermion 
fields [26-29] 

^(±) = P{±)^, 

P(±) = ^(l±7o73) (2.9) 



Since U'^ = 1 the equations (|2.8|) may be rewritten in the form 

^f/"D^^(±) = M(y2?7±7o^(T) - ^.7,*(±)) (2.10) 

From these equations it follows that both fields \E'(-|-) are dynamical fields because both of 
them contain the derivative with respect to the "time" coordinate x^. This distinguishes 
heavy- fermion QED from ordinary QED where the fields ^'(-i-) are not independent [26-29] in 
light-cone quantization: The field \E'(+) can be related to the dynamical field \E'(_) in ordinary 
QED. (See, however, the discussion of zero-modes in [pO| .) 

The Lagrangian (|2.1|) is gauge invariant. This means that the classical theory contains 
"first-class" constraints and we need a quantization prescription for systems with constraints 



such as that provided by Dirac [21 



According to this procedure "primary constraints" arise when one cannot relate a ve- 
locity to its corresponding canonical momentum. The consistency condition of a primary 
constraint, which means that the constraint must be conserved in time, is used either as a 
condition which defines a Lagrange multiplier function or as a "secondary" constraint. There 
could even be "tertiary constraints", and so on. This procedure terminates when no new 
constraints appear. 

All constraints are divided further into two groups. Constraints whose Poisson brackets 
with all other constraints vanish on the constraint surface are called "first-class" constraints. 
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Otherwise, they are called "second-class" constraints. If a constraint is first class, a sub- 
sidiary condition (gauge condition) must be imposed in order to determine the corresponding 
Lagrange multiplier function. When all constraints finally become second class, we can invert 
the matrix of Poisson brackets between the constraints and replace the Poisson bracket by 
the Dirac bracket. The quantization procedure consists of replacement of the Dirac bracket 
by a commutator for bosons or an anticommutator for fermions. 

Let us consider the quantization of the theory in the restricted region 
—L < x~ < L, — L < < L and impose periodic boundary condition for the bosonic 
variables and antiperiodic boundary condition for the fermionic ones. In such a restricted 
region, the fact that the classical equations are first order in the "time" leads to the 
well known problem of zero modes. To separate the zero mode contributions, it is useful to 
generalize the orthogonal projection operators found in [O] 



^(a)*0U) = / V(a){x,y)(l){y)dy, 

V^a) = {V,V„Q,Q,), J = 1,2, 
Q = I-V,Qj=I-Vj, 

V{x,y) = ^S{x±-y^), 

Vj{x,y) = ^ej,6{x'-y'')6{x--y-), (2.11) 

where x = {x~ , x^, x^) = (x~, x±), antidiag ejk = (1, 1), and X is a unit operator: 
I*(j)[x) = (j){x). Here V, Vj are the projection operators into the zero- mode sectors of x~ ,x\ 
respectively, and Q, Qj are the projection operators which eliminate the zero modes in the 
sectors of x~,x^, respectively. We shall speak of V as the projector onto the P-sector, and 
Q as the projector onto the Q-sector. 

In the same way it is possible to define the product of two projectors 

^(a)P(b) * <P{X) = V^a) * (P(fe) * Mx) = ^(6)^(a) * <l){x) (2.12) 

Due to this definition, the operator V1V2 is orthogonal to Qi, Q2 and QiQ2- One can 
verify that the operators V1V2 and Qi + Q2 — Q1Q2 are projection operators. Further, 
^17^2 + Qi + Q2 — Q1Q2 = 1, so that the two subspaces defined by these projection operators 
span the space. 

We define the canonical momenta 11^ as 

_drC 

= d+<p (2.13) 
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Then 



n+(= n^, 



= 0, 

= i_ - d-A+ 







(2.14) 



The velocity can be expressed through the momentum n_ and we have five primary 
constraints ^^'^^ ~ 0, where 



n+ 

Uj - d-Aj + djA. 
- i^JjJU_ 



(2.15) 



The canonical Hamiltonian density on the constraint surface ( |2.15| ) is 
He = -Il^Ul\U+d^ - Ujdj)^ - iMU^I/lUZ^^ 

-ieU^Ul'U'^A^^ + - (a_n_ + djUj)A+ + ^F.^F, 



jk 



Consider the Poisson brackets between the constraints (f)^^^ from ( 2.15 ) 



(2.16) 



(2.17) 



From this expression we see that the constraints (j)^^^ commute only in the P-sector. Hence 

in the P-sector the consistency conditions for the constraints (f)^^^ may lead to secondary 

constraints. Thus let us consider the decomposition of the primary constraints 0+'' and (p^j^^ 
into the Q and V sector constraints 



Q^iUj-d^Aj + djA.), 

p*n+, 

V*(Ilj+djA_) 

Following the Dirac prescription we construct the primary Hamiltonian density 
+Af VJ.^'S) + A^Vr^ + Af 



(2.18) 
(2.19) 



(2.20) 
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where A\. ' , A^- ' are Grassmann-even Lagrange multipliers for the constraints (|2.18|) , 
( 2.19| ) and A,i,,A-^ are Grassmann-odd Lagrange multipliers for the fermionic constraints 
(f)^^\(f)^\ respectively. The consistency conditions for the constraints (f)^^\(f)^^ do not give 
new constraints, and hence lead to the determination of their Lagrange multipliers. We will 
consider the Q and V sectors separately. 



3 Q-sector 



In this Section we are going to consider the quantization procedure in the Q-sector 
In this sector the primary Hamiltonian density is 



(3.1) 



The consistency conditions for the constraints 
and for the constraint 0+'^^ 



J determine their Lagrange multipliers, 
gives a new constraint 



(3.2) 



Q * (a_n_ + djUj + ieHvp*) 

This constraint does not commute with the fermionic constraints. Thus let us consider a 
new "Gaussian law" constraint which is equivalent to ( p.2| ) in the presence of ( p.l5|) 



Q * (9_n_ + djUj + ieU^^ + ie^%) 



(3.3) 



and, on the constraint surface, commutes with all constraints. The fermionic constraints are 
second class constraints and commute now with all nonfermionic ones. Therefore we will 
omit them in future consideration. 

The consistency condition of the constraint ( |3.3|) gives neither a new constraint nor 



determines any Lagrange multiplier. The constraints 



are second class, 



and the (/>+'^'', 0+'^'* are first class. We need two gauge conditions. We choose them in the 
form 



Q * d^A_ 

So in the Q-sector the constraints $(2) are 



(3.4) 



$(S) 



-Q * Ali(9_n_ + dkUk + ze(U^^ + ^%)) 



4^) 



(3.5) 
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where A_ = and is a nondegenerate operator whose matrix elements in the Q-sector 



are 

(A:'),,j = H^[x- -v-)i{,x^-Vi), 

H-(x-y) = ^_(£__^ (3.6) 

The constraints $(2) are second class. The matrix of Poisson brackets between these 
constraints 

Aa,b{x,y) = {(i)i^\x),(t)f\y)}, a,b = l,...,6, (3.7) 

is quasi-diagonal (block-diagonal) and its non-vanishing elements are found to be 

Ai,2U-/y) = -A2,i{x,y) = -Q{x,y), 
A3,4{x,y) ^ -A4^six,y) = -Qix,y), 

Aj+4j>+4{x,y) = -25jfd-Q{x,y) (3.8) 

The inverse matrix A~^{y,z) 

'dyA,4x,y)A-l{y,z) = 5,,6^Q\x-z), (3.9) 
6^^\x -z)^ D^{x- - z-)S{x^ - z^), 

dQ{x-z)^S{x-z)-^, (3.10) 

has non- vanishing elements 

^rjfey) = -^2,1 fe^) = -Q{x,y), 
All{x,y) = -Ajj(x,y) = -Q{x,y), 

^7+4j'+4fey) = -2^jj'^i^± - z±)G^{x'- - z~) (3.11) 
G^{x-z) = '-^^-^ (3.12) 

The functions D^{x — z) and G^{x — z) are the delta function and the matrix element of 
the operator dz^, respectively, in the Q sector. 
Consider two variables 

u;f^ ^Q*Aj (3.13) 

It is easy to show that the Dirac brackets between these variables are 

{ujf\x),u;f\y}n = -^<5,y5(:r^ - y^)G'^ix- - y-) (3.14) 
Thus we can consider the variables o;]®^ (x) as the physical variables in the Q-sector. 
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4 P-sector 



In this Section we will consider the canonical quantization procedure in the P-sector. 
The primary Hamiltonian density in this sector is 

= K + ArvS'"^ + y^fUr^ (4.1) 

The P-sector is a space if^'^) of functions which depend on the variables X"*" only. On 
the other hand the P-sector is the sector of zero-modes in the ^"-direction. To take into 
account the zero- mode contributions in the directions x^^x^ one can decompose the space 
into a direct sum of two orthogonal subspaces 

lj{V) ^ jj{VViV2) jj{V{Qi + Q2-QiQ2)) (^4 2) 

where if (^^1^2) is the space of zero-modes in all directions and is defined by the projector 
PP1P2 and the space if(^(2i+S2-2i22)) is defined by the projector P(Qi + Q2 - QiQ2)- 
Such a decomposition leads to the corresponding decomposition of the primary constraints 
4>^^''^\(f)f''^^ ( p. 191) , and the primary Hamiltonian density Ti'-^''^-' (|4.1| ) 

<t/r = (^f o . (4.3) 

= Ho + Ai;VV" + Aj'ri" + a5'V« + AfV'f (4.4) 

Here A™, A?,Af' ,A<" are the Lagrange multipliers; the primary constraints in the PP1P2- 
subsector are 

= PP1P2 * n+, 

= pp,P2 * n,- (4.5) 



and 



V'+^=P(Qi + Q2-QiQ2)*n+, 

= P(Qi + Q2 - Q1Q2) * (n, + djA_), (4.6) 

are the primary constraints in the P(Qi + Q2 — QiQ2)-subsector. 

The consistency condition for the constraints (p^^\(p^^^ and do not determine 

their Lagrange multipliers, and instead give new constraints in the form (we denote them as 
X+\xf\x+\xf\ correspondingly ) 

= PP1P2 * (zen*^), (4.7) 

xf'> = VV1V2 * i-ieU^UZ^Uj'^) (4.8) 

= P(Qi + Q2 - Q1Q2) * {d,U, + leU^^), (4.9) 

xf = ViQi + Q2- Q1Q2) * (5,n„ + dkFkj - zen^t/li?7,^) (4.10) 
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As it happened in the Q-sector the constraints ( |4.7| )-( |4rT0| ) do not commute with the 
fermionic constraints. Thus let us consider new constraints which are equivalent to (|4.71 )- 
( |4.1CI| ) in the presence of ( |2.15D 

X^l^ =VViV2*p, (4.11) 

xf =VViV2*Pj. (4.12) 

= ViQi + Q2- Q1Q2) * (d.Uj + p), (4.13) 

xf = V{Qi + Q2 - Q1Q2) * (5,n_ + dkFk, - p,), (4.14) 

where 

p = ze(nvp^+^%), (4.15) 

and 

Pj = ^P, (4.16) 

which is the consequence of the heavy mass limit. 

Let us consider the quantization procedure where the strong constraint 

VV1P2 *p = (4.17) 

holds. The quantity VV1P2 * p is the total light-cone electric charge of the fermions per unit 
volume. This is proportional to the usual electric charge, which must vanish in a compact 
system as a consequence of Gauss's law. 

The consistency condition for the constraints ( |4.13| ),( ^.14| ) does not give any new con- 
straints. So in the P-sector we have fourteen constraints : (p^^ , 99^"^'', x+\ xf\ '^+\ i^f\ 

X^+ 1 xf \ 0^'*- The constraints p)l\ ipf\ and are first class and the rest are 
second class. We need four gauge conditions. 

We consider the gauge conditions for the first class constraints ^9+'* and ip^^^ in the form 

<^G = VV1V2 * A+, 

tpS^ =V{Qi + Q2-QiQ2)*A+ (4.18) 

Let us discuss the gauge conditions for the first class constraint in the V{Qi + Q2 — 
QiQ2)-subsector. Note that in this subsector the operator 

A = dkdk (4.19) 

is an invertible operator. Let us choose the gauge condition for the constraints x+ the 
form 

x'g^ = V{Qi + Q2 - Q1Q2) * ^-'duAu (4.20) 
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The constraints xp do not commute with the primary constraints ip^p but do commute 
with the rest. On the other hand, the primary constraints ifj^^ do not commute with the 
gauge condition x^g - Thus let us consider a hnear combination of the constraints ifjf''^ and 

(2) 

x+ 

) = i^f^ + a,xf (4.21) 

in such a way that ip'-^^ will commute with the gauge condition x^q ■ Here aj are some 
functions to be found. We obtain 

= p(Q, + Q2-QlQ2)*(n,^-|p + fi',-4-), (4.22) 

n,^ = (^..-^)n. (4.23) 

Here Hj- are two-dimensional transverse momenta. Let us introduce longitudinal momenta 
and an analogous decomposition for the potentials 

n, = h^^ + hJ, hI! = ^h,, (4.24) 

Aj = Aj -\- A^j , 

Af = iS,,~^)A,, A^] = ^A, (4.25) 

With this decomposition the constraints depend only on transverse or longitudinal com- 
ponents, but not both 

X? = P(Qi + Q2- Q1Q2) * {d,u] + p), (4.26) 
Xg =V{Q, + Q2- Q1Q2) * A-'dkAl (4.27) 
xf =V{Qi + Q2- Q1Q2) * (c',H_ + ^Aj - pj) (4.28) 

Consider the primary first class constraints <^f^- Naively one might try to choose the 
gauge condition as VV1P2 * Aj ^ 0, which means that the zero-modes of Aj are eliminated. 
This contradicts, however, the fact that the integral over a closed loop 

Ajdxj (4.29) 



is gauge invariant |]3T|, and need not vanish if taken over a non-contractible loop. Instead, 
we choose the gauge condition Lp^jQ ^ in the form 

^%=VV,V2*A,-f,{x^), (4.30) 
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where fj{x^) is any function of the "time" x'^ . The constraints (p^^Q depend on "time" 
exphcitly. It is well known that the consistency condition for this type of constraint can be 
written in the form 



(4.31) 



This relation determines the Lagrange multiplier A^^'' in terms of the derivative d+fj. In the 
next Section we will show that the physical Hamiltonian does not depend on the function 

So the constraints <I>'^) ss in tlie V sector are (we omit here the fermionic constraints) 



i-i> 

(2) 

X g -- 

(2) 



4? 



> PPiP2-subsector 



(4.32) 



>V{Qi + Q2- QiQ2)-subsector 



The constraints ( [4.32|) are second-class constraints and the matrix of their Poisson brackets 
is block-diagonal. 

The constraints 4>'^a\ (a = 1, ■ ■ ■ , 10), determine the nonphysical variables. The physical 
variables are among the constraints 0j+io, 0}+i2- The matrix of Poisson brackets between 
these constraints is antidiagonal and its nonzero elements are 



{0i?iofe),0r+12(y)} 



{V) 



A 



j+10J'+12 



= -6j,,A{V{Qi + Q2- QiQ2))ix,y), 
A+i2,i'+io(2I,y) = 5,yA(P(Qi + Q2-QiQ2))(x,y) 

The inverse matrix is antidiagonal too, and it easily found to be 



(4.33) 
(4.34) 



^i+i2j'+io(2ii y) 



^i+io,i'+i2fe y) 



-5n'^inQi + Q2-QiQ2)){x,y), 



QiQ2)){x,y) 



(4.35) 
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Consider two pairs of conjugated variables 

JP = V{Qi + Q2 - Q1Q2) * n_ (4.36) 
=VViV2*A^, 

J.P ^ VV1V2 * n_ (4.37) 
The non-zero Dirac brackets between these variables are 
{ujPix),jPiy)}n = 2{V{Qi + Q2-QiQ2)){x,y) 

= ^(5(xx-y^)-^), (4.38) 

{u:f\x),u:f\y)}n = {VV,V2){x,y) = (4.39) 
Thus one can consider these variables as physical variables in the P-sector. 

5 Physical Hamiltonian 

We have come to the step where we can find the "physical" Hamiltonian. Such a 
Hamiltonian is defined as 

Hf^ = Kk.o (5.1) 

where 

$ = ($(S), $(^), (5.2) 

are all constraints in the problem under consideration. 

From eq. ( |5.1| ) it follows that to construct the Hamiltonian density T-C^^^ we have to express 
all variables through the physical variables using the constraints $ «i 0. For this purpose 
let us return to the expression ( p^.l6 ) for the canonical Hamiltonian and rewrite it in the 



following form 

H^ = J nPdx + J n'.dx (5.3) 

Here 

nP = '^JJI{-iU+d^ + iUkdk)^ + M^^ (5.4) 
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is the Hamiltonian density of fermions. The expression Ti'^ 

K = \{Q * n_)2 + ^{vViV2 * n_ + viQi + Q2- Q1Q2) * n,)^ 
+(Q * n^)ivViV2 * n_ + viQi + Q2- Q1Q2) * n_) - 

(9_n_ + djUj){Q * A+ + VV1V2 * A+ + ViQi + Q2- Q1Q2) * A+) 
+\iQ * Fjkf + \{V{Qi + Q2- Q1Q2) * F,kf + 

+e^^^^(Q * + VV1V2 *A^ + V{Qi + Q2- Q1Q2) * A^) (5.5) 

is the Hamihonian density of the electromagnetic field and its interaction with the fermions. 
Then the physical Hamiltonian becomes 

-Hf^y = Ki^^o (5.6) 

Due to the property 

e^J/IU^^ = U^p (5.7) 

and the strong condition ([4.17|) one can neglect the term e'^Jj}U^'^VViP2 * A^ in the Hamil- 
tonian density H'f^y' (U). This is because the physical Hamiltonian H''^^^^ , which corre- 
sponds to the density lUf'^'^^, by itself is proportional to the zero-mode of Ti'f'^'^^ 

Thus the physical Hamiltonian does not depend on the functions fj{x^) involved in the 
gauge condition ( [4.30| ). 

Using the constraints $ ^ we are able to write 

Q * H_|*^o = -{djoof^ + eQ* dZ^^W-^). 
Q * A±ksiO = 0, 

VViV2*A+\^^o = Q, 

ViQi + Q2- Q1Q2) * ^+1*^0 = 0, 

VV1V2 * Aj\^^o = fj{x+), 

V{Qi + Q2 - Q1Q2) * 4.1*^0 = 0, 

P(Qi + Q2- Q1Q2) * 4k.o = + 

V{Qi + Q2- Q1Q2) * Fjkl^^o = V{Qi + Q2- Q1Q2) * iffek^o 

= Uz\Ukd,-U,dk)JP, (5.9) 
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After substitution of eqs.( |5.9| ) into eq.( |5.6| ) we get 



UUr, 



(S) 

jk 



(S) 

jk 



jk 



(P) 



where 



(S) 



jk 



(V) 



dj^k 

T-l 



(S) 



(5.10) 



dkUJ 



(S) 



jk 



(V) 



(5.11) 



Having the physical Hamiltonian and the Poisson brackets between all physical variables 
one can now quantize heavy fermion QED. According to the Dirac prescription the quanti- 
zation procedure consists of replacing the Dirac brackets of bosons with bosons (or bosons 
with fermions) by a commutator and that of fermions with fermions by an anticommutator 



{A, B}£i —i[A, i?]-for bosons, 
{A,B}o —i[A, B]^-foY fermions 



(5.12) 



where the symbol ' means that the classical object A becomes an operator A. Then we 
come to the following equal "time" commutation relations 



{V) 



i 



--5jk5{x±_ - y±)G{x - y' 
4(5(^x-yx)-^), 



\UJ 



8L3' 



^5ix - y) 



(5.13) 
(5.14) 
(5.15) 
(5.16) 



LJa, where 



The physical Hamiltonian H^^y^ can be obtained from ( |5.1CI| ) by the change uo, 
ijja are all operators satisfying the relations ( ^.13|) -( pn6D . 

We wish to find a realization of the commutation relations ( 5.13| )-( ^.16 ). Using the 
periodic boundary conditions for and eq. (|5.13|) one may write (from now on we will 

omit the operator symbol " ) 

1 ^ 1 



J^\x) 



E 



2PJ 



(aj[P_n) exp(-iP„ ■ x) + a+(P„) exp(iP„ ■ x) 



(5.17) 
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where 



= n+ = 1,2,..., n^' = 0,±1,±2,..., Q = 8L^ 



(5.18) 



and for any vectors a-h = a^h —ajbj. The operators aj{P_n) and aj{P_n) satisfy the following 
nonzero commutation relations 



'jj'Op p 
J J —n—17 



(5.19) 



Consider the operators (^i^^(xx). Their decomposition into a Fourier series can be written 
in the following form 



'^r^(a;±) = -4= exp(-iP_j_„ ■ x±) + a+(Pj_„) exp(jP_j_„ • x± 



1 



)) 



(a-a+). 



V" p,„ 



(a + a+). 



where 



TT 



-P±n = -^"-^ rv' = 0, ±1, ±2, . . . , a_L • &_L = Oj&j, 



(5.20) 



(5.21) 



and the operators a(Pj_n), a"'"(P_Ln), a, a'^ satisfy the following nonzero commutation relations 



a(P±„),a+(P_j_„) 
[a, a"*"] = 1 



(5.22) 



(v) 

The realization for the zero-mode operators (^3 4 is 



<^3 



(p) 



;c-c+). 



[c, c+] = 1 



(5.23) 



The operators aj(P„), a(Pi„), a, c and a^(P„), a"'"(Px„), a+, c"^ can be interpreted as annihi- 
lation and creation operators. 

The same decomposition in terms of creation and annihilation operators can be made for 
the fermion operators ^'(^), ^'(^) satisfying the antiperiodic boundary conditions 

*(x) = Yl ihiQn)u+ eM-iQn ■ x) + d^{Qn)ut exp(iQ^ • x)) (5.24) 
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where 

Qn = iQt, Qi) = f (^^ + 1/2, + 1/2), 
n+ = 0,1,2, n^' = 0,±1,±2,--- 
A±4 = 4, A± = ^(l±^, 

= ±(5„^, = 0, (5.25) 

and the operators hp{Q^,d'p{Q^ satisfy the following nonzero commutation relations 

MQJX{QJ]+ = mQJ,dl{QJ]+ = 6.p6q^q^ (5.26) 

6 Calculation of Heavy Fermion Potential 

We wish to see the physical Hamiltonian H^^'^'^ found in the previous Section in action. 
The example we choose is the static heavy fermion potential. (This calculation is the analog 
in QED of the calculation of the quark-antiquark potential in QCD.) For this purpose we 
consider old-fashioned perturbation theory up to second order 

E ^ E(o) + eE(i) + e2E(2) (6.1) 

The physical Hamiltonian H^^^^ can be expressed in terms of the normal product of the 
creation and annihilation operators introduced in Section 5 

HP^y^operator) = H(^o) + ei/(i) + e^H(^2) (6.2) 

The construction of the normalized eigenstate which describes the two static heavy 
fermions r apart yields 

M M 

Iff) = l)-'/' E exp(-^(g„ - MU) ■ r)ht{V2M - Q+, -Qxn)4(QJ|0) (6.3) 

Q Q 

— m — n 

where a ^ (3 and we choose 

[/+ = [/- = [/I = [/2 = 
v2 



In eq. ( |6.3| ) we have the effective summation over which has modulus less than M. It 
can be shown that the corresponding eigenvalue (-^{o)) of H(q) turns out to be \/2M and the 
mean value of the operator H{i) in the state ( |6.3|) is zero 

E(i) = (/7|i/(i)|//) = (6.4) 
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The last term in eq.( |6.1| ) is 



(6.5) 



where |m) is an eigenstate of if(o) and E"^-^ is the corresponding eigenvalue and the summation 
does not cover the state \ff). After some straightforward calculation we get 

(/7|i^(2)|//) = E + cos(P„, ■ r)] (6.6) 



and 



^\{ff\HnM\' 1 ^ 1 P_ 



2 



m E(0) - Ji(o) ^ + P±n Pn 

If we remember that r = f"^) and make use of the new vectors 

P„ ^ (P„\ P„^ v^P+), r ^ (r\ r^, r,^) (6.8) 
then (in the limit L — > cxd), 

^ 1 _ 



»=^<^. = S pill + ■ = - j;!^ + Coust. (6,9) 



where is the Cartesian z-component of the vector r. The additional factor l/\/2 in the 
expression ( |6.9|) comes from the fact that the energy e^E(2) corresponds to the Hamiltonian 
P~ in the light-cone coordinates. In the reference frame where P = we have 

P- = ^po 

This gives the potential in the Cartesian coordinates 

Vir) = -— (6.10) 

and clearly shows the Coulomb potential together with an irrelevant infinite constant in 
eq.( |6.9|) (which is proportional to M) arising from the fermion self energy. 



18 



7 Conclusions 



We have considered the hght-cone canonical quantization of the heavy fermion QED 
taking into account the zero-mode contributions exphcitly. We have imposed periodic bound- 
ary conditions for bosonic fields and antiperiodic ones for fermionic fields. As in ordinary 
QED, this model is gauge invariant, which means that there are unphysical degrees of free- 
dom. To quantize the theory, we have used the Dirac algorithm for a system with first- and 
second-class constraints and the corresponding gauge conditions. In order to make explicit 
the role of the zero-modes, we have considered gauge fixing and quantization procedures in 
the zero-mode and non- zero-mode sectors, separately. In all sectors we obtained the physical 
variables and their canonical (anti-)commutation relations. The physical Hamiltonian was 
constructed by excluding unphysical degrees of freedom. We have considered the role of all 
physical fields in the calculation of the potential between static heavy fermions. 

We suggest that this approach can be used for the case of finite-mass QED or, perhaps, 
one might first address heavy quark QCD. An important goal is the calculation of the QCD 
quark-antiquark potential. 
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